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Surface plasmons (SP) in a semi-bounded quantum plasma with degenerate electrons 
(e.g., a metal) is considered, and some interesting consequences of electron Pauli block- 
ing for the SP dispersion and temporal attenuation are discussed. In particular, it is 
demonstrated that a semi-bounded degenerate plasma with a sharp boundary supports 
two types of SP with distinct frequencies and qualitatively different tempora l attenua- 
tion, in contrast to a non-degenerate plasma that only supports one type of SP (| Guernsey 
1963). 



1. Introduction 

The properties of surface plasmons (SP) in bounded plasma structures are defined, 
among other things, by the dielectric properties of the plasma medium that sustains 
them. The latter are often (e.g., in metals, for which the electrons are strongly de- 
generate) significantly affected by the quantum nature of the charge carriers and their 
interaction in the medium. This can affect the properties of SP in a non-trivial way, via 
modification of analytic properties of the medium response. In particular, quantum ef- 
fects (du e to Pauli blocking and o verlapping wave functions of free charge carriers in the 
medium (IShukla fc Eliasson|[201ll)), when significant, can n iodify the dispers ion, damp- 
ing (jTyshetskiv et aLll2012M ) and spatial attenuation of SP (| Vladimirovll 1 9941 ) supported 
by a bounded quantum plasma. 

In this paper, we show another consequence of quantum degeneracy of plasma electrons 
on SP properties, exemplified by a simple case of SP in a semi-bounded collisionless 
plasma with degenerate electrons. In particular, we show that such system supports two 
types of SP, with different frequencies and qualitatively different temporal attenuation, 
in cont rast to a case of non-degenerate semi-bounded plasma that only supports one type 
of SP (|Guernsevlll969h . 



2. Method 

2.1. Model and assumptions 

We consider a semi-bounded, nonrelativistic collisionless plasma with degenerate mo- 
bile electrons (Te <^ ep, where Te is the electron temperature in energy units, = 
h^{2>iT^nef'^^ /2me is the electron Fermi energy), and immobile ions; the equilibrium 
number densities of electrons and ions are equal, noe = Jioi = (quasineutrality) . The 
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plasma is assumed to be confined to a region a; < 0, with mirror reflection of plasma 
particles at the boundary x — separating the plasma from a vacuum at a; > 0. 

We will look at SPs in the non- retarded limit, when their phase velocity is small 
compared with the speed of light. In this limit, the SP field is purely electrostatic, hence 
we can restrain ourselves to considering o nly electrostatic oscillatio ns in the considered 
system. Following the discussion of Ref. ( Tvshetskiv et a/] 12012 we adopt here the 
quasiclassical kinetic description of pla sma electrons in terms of t he 1-p article distribution 
function /(r,v,t) = /(x,r||,W2:,V||,t) (jVladimirov fc Tvshetskivll201l[) (where r|| and V|| 
are, respectively, the components of r and v parallel to the boundary, and x and Vx are 
the components of r and v perpendicular to the boundary), whose evolution is described 
by the kinetic equation 

df df e fdcjy df ^ dcj, df\ „ ^2 1) 



di 

dt 



ox " 

where the electrostatic potential 



9r 



uc,ri 



Aire 



^dx dvx dr\\ Svy ' 
, t) is defined by the Poisson's equation 

/(r,v,i)d3v-no 



(2.2) 



In the absence of fields, the equilibrium distribution function of plasma electrons /o(v) 
is defined by an isotropic Fermi-Dirac distribution, which in the limit <^ ep (where 
is the electron temperature in energy units, and ep — {Sn^hPrie)^^^ /2me is the electron 
Fermi energy) reduces to 



/o(^;) 



2ml 
(27rft)3 



1 + exp 



2m?' 

^^^(7[z;i.(ne)-^;], (2.3) 



where vp(n^ ~ ^J2ep{ne)/me is the electron Fermi velocity, <t{x) is the Heaviside step 
function. 

The condition of mirror reflection of plasma electrons off the boundary at x = implies 
fix = 0,r||,-U;j., V||,t) = f{x = 0,r||,w^, V||,t). (2.4) 

2.2. Initial value problem 

We now introduce a small initial perturbation fp{x,Y\\,Vx,'v\\,t = 0) to the equilib- 
rium electron distribution function /o(w), |/p(a;, r|| , ti^;, V|| , i — 0)| <C fo{v), and use 
the kinetic equation (|2.ip to study the resulting evolution of the system's charge den- 
sity p(x,r||,t) = e [/ /(a;, r|| , V, t)(i'^v — no] , and hence of the electrostatic potential 
4>{x, r|| , t) defined by (j2.2l) . Introducing the dimensionless variables — uj/ujp, K — kA^, 
V — v/vp, X = x/Xp, R || = r ||/AF, A_f = vp/\/3ujp, ujp — (47re^no/we)"^/^, and follow- 
ing Guernsey ( Guernsevlll969l ). the solution of the formulated initial value problem for 
p(X, R||,T) with the boundary condition (|2.4|) is 

p(X,R||,T) = eno/5(X,R||,r), (2.5) 

where 

1 



p(X,R||,T) 



Pk(T) = ^ / 



(27r)3 



dK^ 



d^Kn 



Pk(r), 



p{n,K) e-''" dn, with cr > 0. 



(2.6) 



(2.7) 



The integration in (|2.7|) is performed in complex n plane along the horizontal contour 
that lies in the upper half-plane Ini(r2) — a > above all singularities of the function 
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p{^l, K). The function p{il, K), defined as the Laplace transform of pk{T): 

/•oo 

Jo 



(2.8) 



is found to be 



e{VL,K) 



(PY- 
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G(V,K) 
V3K- V 
1 



e{n,K) 



dKL 



.3v^iX^^2.9) 



where the Fourier transforms G(V, K) and G'(V,K') of the (dimensionless) initial per- 
turbation are defined by 



with 



/ + 00 
-OO 

g(X,V;,V||,K||) - yd^Rii e-*K|,R„ /^(x, R,, , V^, V,, , 0), 

3 

/p(X,R||,y,,V||,0) = ^/p(X,R||,F,,V||,0), 

II II 

where iiTii = \K\\\, K ^ \K\, K (X^,K||),iC' = |K'|, and K' = {K'^,K\\ 
e{n,K) and C(^^,^i^||) in are defined (for Im(f]) > 0) as follows: 



1 



K.a/o(v)/av^3 



V^K^ J - \/3K • V 



2tt 



+ 00 



(2.10) 

(2.11) 

(2.12) 
The functions 

(2.13) 
(2.14) 



with 



/o(V) = -^/o(V) = -^/o(v)U^„^v- 

uq no 



For fully degener ate plasma with electron distributio n (12. 3p . the function £(51, K) becomes 
(for Im(n) > 0) (|Gol'dmanl[l947l : lAlexandrov et fl/.lir984D : 



1 - 



In 



, lm{n) > 0, 



(2.15) 



where ln(z) is the principal branch of the complex natural logarithm function. 

Note that the solution p.9p differs from the corresponding solution of the transformed 
Vlasov-Poisson system for infinite (unbounded) uniform plasma only in the second term 
involving (^{fl, K\\)] indeed, this term appears due to the boundary at x = 0. 

The definition (|2.8p of the function p{n, K) of complex ft has a sense (i.e., the integral 
in (|2.8[) converges) only for Im(r2) > 0. Yet the long-time evolution of Pk{T) is obtained 
from (|2.7p by displacing the contour of integration in com plex fl plane f rom the upper 
half-plane Im(il) > into the lower half-plane Im(r2) ^ (|Landaulll946f ). This requires 
the definition of p(ri, K) to be extended to the lower half-plane, Im(n) ^ 0, by analytic 
continuation of (12. 9p from Im(f7) > to Im(ri) ^ 0. Hence, the functions 



J(fl, K) EE / d^V 



GjV.K) 

n - V3K ■ V ' 



(2.16) 
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Figure 1. A sketch of the deformed integration path in (12. 7p in complex Q plane, with contri- 
butions of poles (solid circles) and branch cuts (dashed lines) of 1/({Q,, K^^). The singularities 
due to l/e{Q, K) are not shown, but they also contribute to (|2.7p . yielding volume plasmons. 



e{^l,K), and ({n,K\^) that make up the function p(f2,K), must also be analytically 
continued into the lower half-plane of complex fl, thus extending their definition to the 
whole complex D, plane. With thus continued functions, the contributions to the inverse 
Laplace transform (|2.7p are of three sources ( Guernsey 1969[ ): 

(a) Contributions from the singularities of /(Jl, K) in the lower half of complex ft plane 
(defined solely by the initial pe rturbation ^(V , K)); with some simplifying assumptions 
about the initial perturbation ( Guernsey 19691 ) these contributions are damped in a few 
plasma periods and can be ignored. 

(&) Contribution of singularities of 1/£{^1,K) in the lower half of complex plane, 
of two ty pes: (i) r esidues at the poles of l/e{^l, K), which give the volume plasma os- 
cillations (| Guernse y 1969), and (ii) integrals along branch cuts (if any) of l/e{n,K) in 
the lower half-plane of com plex f^, which can lead to non-exponential at tenuation of the 
volume plasma oscillations (|Hudson 1962 : Krivitskii fc Vladimirovl [19911) . 

(c) Contribution into (|2.7p of singularities of 1/^(17, K\^) in the lower half of complex fl 
plane, of two types: (i) residues at the poles of K^^), correspondii ig to the surface 

wave solutions of the in itial value problem in the considered system ( Guernsevlll969l : 
Tvshetskiv et al. 201261 ). and (ii) integrals along branch cuts (if any) of 1/C(r2, isTii) in 



the lower half-plane of complex VL. 

Below we consider the latter contributions from poles and branch cuts of 1/C,{Q.^K\\) 
in the lower half-plane of complex fi, as illustrated in Fig. [l] and show that they yield 
two types of electrostatic surface oscillations with different frequencies and qualitatively 
different temporal attenuation. 



3. Two types of surface oscillations 

3.1. Contribution of poles of \/C,{^,K\\) 

The contrib ution of poles of IJCi^, K \\] into (|2.7p leads to exponentially damped surface 
oscillations ( Tvshetskiv et al. 201261 ) 

p(P°"''')(T)cxe-|r=l^cos(f7,r), (3.1) 

with frequency = ujg/ujp and damping rate T s = ^s/^p obtained from the dispersion 
equation C,{yi,K\\) — 0. The frequency asymptotes are 



1 + 0.95X||) , for K|| < 1 



V2 

V3K« f 1 + 2 exp 



-2 - iKf, 



for Kii ^ 1 (zero sound) 



(3.2) 
(3.3) 
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3 . Integrate the 
analytically continued 
approximate p(fl,K) 
along the branch cut, 
^ near the branching point y 



1 . Approximate p(n,K) in the upper 
vicinity of the branching point flX^w)^ 
via elementary functions of H — Qy 
I (namely, ^JQ — Q,, - a 2-branch function) 



Re(n) 



2. Analytically continue the approximate 

p (Q, K) fi'om the upper to the lower 
vicinity of r2,,(/i'||), choosing the cut of 



- Q„ as shown 



Figure 2. Evaluation scheme for the integrals along branch cuts, as outlined in Sec. 13.21 



The absolute value of the damping rate is a nonmonotonic function of K\\. At small K\\, 
it increases linearly with , 

|r,(if||)| «2.1\/3-10-2if||, (3.4) 

reaches maximum jr^l « 6.2 • 10"'^ at sa 0.4, and then quickly decreases at > 0.4. 
Since the maximum growth rate is small, the su rface oscillations due to the poles of 
l/({fl,K^\) are weakly damped at all wavelengths (jTyshetskiv et a l. 20126,). 

3.2. Contribution of branch cuts of 1/(1{Q,,K\\) 

For degenerate plasma, the analytically continued function C,{Q,^K\\) has two branching 
points on the real axis of the complex Q, plane ai Vt — ±Q,y{K\\), where Vl^{K\\) G R is 
the solution of equation 

£(l],i^ll) = e(r!,i^)|^^^o=0, (3.5) 

with the corresponding branch cuts going down into the Im(i7) < part of the complex J7 
plane, as schematically shown in Fig. [1] (see Appendix lA)) . Let us consider the contribution 
of the integration along these branch cuts into the inverse Laplace transform (12. 7p . The 
branching points lie above the poles fi^ — i\^s\ of 1/C (since the latter lie below the real 
axis of the f2 plane), therefore we can expect the contribution of the integration along 
the branch cuts into (|2.7p to be at least as important as the contribution of the poles, if 
not to exceed it. 

At large times T ^ 1, the main contribution into the integrals along the branch cuts 
comes from the small vicinity of the branching points, so it sufhces to approximate the 
second term of (|2.9p near the branching points in the lower semiplane of complex Vl. This 
can be done in two steps: 

(a) approximate /5(r2,K) defined by p.9|) in the upper vicinities of the branching 
points, in terms of elementary functions; the approximate function should have the same 
branching points as the original one; 

(6) analytically continue these approximations into the lower vicinities of the branching 
points, choosing the branch cuts to go down. 

Then we can perform the integration of thus obtained approximations along the branch 
cuts in the vicinity of the branching points. This scheme is sketched in Fig. [5] 

The function (j2.9p in the upper vicinity of the right branching point = +i7„(_ftr||) 
can be approximated as (see Appendix [B|) 

/(f^„,K) 



p(fl,K) 



e(Sl,K)\ 



e{n,K)\ 



(3.6) 



1 + ^a{K\\)l3{K\\)K\\y/rr^, 
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where /(O, K) is defined in Eq. (|2.16p and is assumed to vary slowly near the point 
f2 = fit, and to not have branch cuts. The expansion of e(il, K) near = $7^ is 

£(17, K)\^^^^ = e{n,,K) + p{K){n - a„) + o[{n- , (3.7) 

with a{K\^), l3{K\i^) and I3{K) defined in Eqs (|B 3I) - (|B 51) of Appendix |B] (note that in 
general I3{K) ^ (3{K\\), since K 7^ K\\). The approximation (|3.6p with p.7p is expressed 
in terms of elementary functions of J7 — 57^, , which can be analytically continued into the 
lower vicinity Im(r2) < of the branching point — When doing so, the complex 
function VTT^^H^ should be defined so that its branch cut goes vertically down from its 
branching point f2 = 57„. This branch cut can be parametrized as 

ficut = +^v{K\\) -ir], 77 ^ 0. 

The integral around the branch cut in the vicinity of the right branching point f2 = +r2u 
is (the (+) superscript denotes the right branch cut) 

/"OO 

Pk\t) = ^e-^""^^ dr^e-^^[p^+\n,-tf^,K)-p^+\n,-trj,K)\, (3.8) 

where p^^]j,(f2,K) are the left and right branches of the analytic continuation of p.6p 
into the lower semiplane Im(f2) < 0. 

We first assume that the function £{n,,K), analytically continued to Im(r2) ^ 0, does 
not have branching points in some (perhaps small) vicinity of the point fl — fl^; hence 
we have 

near fl — fly. Then the only function in /5^^]j(i7, K) with a cut is \/Q — fly. Choosing its 
cut as specified above, we have 



y/n-fiyj ^ i^^/n - iiyj ^ = ^/\n-ny\e-"/^ = ^\e-'''i^. (3.9) 

Using (1321), (ISH) and (IS21), from we obtain 



.(3.10) 



drj- 



1 



einy,K)-ir]l3iK) + 0{7j^) 



The long-time asymptote of p'^\t) depends on whether £{fly, K) tends to zero or not, 
which in turn depends on the value of K^. Below we consider the two cases: (i) ^ 0, 
so that eifly.K') ^ 0, and (ii) 0, so that e{fly,K') — >• [since fly is the root of 

p.Sp ]. In these cases, p.lOp gives, respectively: 

T3/2 

Similarly for the contribution of the left branch cut, fl = —fly — irj, -q ^ 0, we have 
/'k ^(-^) given by Eqs p.lll) ~ p.l2p with fly replaced with — fly. The total contribution 
of both branch cuts p^"*"^(T) p^'(T) + p^^(T) is then 

~pt''\T)^^-^^^j^ + 0{T-'^'% K^^O, (3.13) 



4+^(T)(x^^;^ + 0(r-5/2), K^^O, (3.11) 
PK\T)^^-r^ + 0{T-'/^), K^~^0. (3.12) 
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(3.14) 



Note that the frequency of these oscillations is equal to the frequency of volume plasma 
waves with the same wavelength, K = K^^ , and thus exceeds the frequency ils of the 
surface oscillations p.ip due to the poles of 1/C. 

Above we have assumed that the function e(f2, K), analytically continued to Im(f2) ^ 0, 
does not branch in at least some vicinity of the branching point fl^ of (^(fl,K\\). How- 
ever, this assumption is violated in a special case considered below. Indeed, the func- 
tion e{il, K) itself, when analytically continued in to Im(rj) ^ 0, has branching po ints at 
r2 — ±y/iK £ M, with the branch cuts going down ( Vladimirov &: Tvshetskivll201lh . Thus 

for Kx ^ ^1{K\\) — the branching points of £(ri, K) merge with the branching 

points of C(^j ^11 )' ^'^^ their respective branch cuts merge at least in some lower vicinity 
of the coinciding branching points. In this case eL{^,K) ^ eR{il,K), and the above 
calculation is modified; instead, we have for Im(il) ^ 



sn{n,K) = l + — 



1 - 



n 



2V^K 

ITT 17 



log 



(3.15) 
(3.16) 



(we still assume that /(r2,K) does not have branching points). Then, after some calcu- 



lation, we obtain for p^^^(i7„ — i?7,K) 

p[+^ {n, - 111, K) - - 111, K) = 



ir], K) in ([3 



3-K 



in this case: 
-HO(mi/2). (3.17) 



Carrying out the integration in p.Sp and adding the similar contribution of the left cut, 
we finally obtain for the contribution of branch cuts in this special case: 



,(cuts)(y)^COs(».r) 



o 



for Kr 



(3.18) 



which has the same frequency as p.l3p - p.l4p . but a different temporal attenuation 
exponent. 



4. Discussion 

We thus see that our system supports two distinct types of surface oscillations, with 
different frequencies and temporal attenuation: 

(a) exponentially damped surface oscillations p.ip with frequency Q,s{K\\), due to the 
poles of p(r2,K) ( Tvshetskiv et al. 2012&); 

[h) power-law attenuated surface oscillations (|3.13p . p.l4p . and p.lSp with frequency 
J7t,(i^ll) > D,s{K\\), due to the branch cuts of p(r2, K). Since the power-law attenuation 
is slower than the exponential attenuation, these oscillations should become dominant 
at large times, and should theoretically become observable. However, by the time this 
happens, the level of the oscillations may well become too low due to the attenuation, 
making a practical detection of this branch of oscillations rather difficult. 
It is interesting to note that, as seen from p.l3p . (I3.14p . and p.lSp . different com- 
ponents in the wave packet making up the surface oscillation of this type is attenuated 
at different rate. Since T~^/^ decays slower than or T~^/'^, the small-ifj, part of the 
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wave packet becomes dominant over the large-Kx part at large times, which corresponds 
to penetration of the perturbation away from the surface and deeper into plasma. 



The presented analysis relies on several assumptions discussed in detail in Ref. (jTvshetskiv et al. 



201261 ) , of which perhaps the most critical ones are the assumptions of collisionless plasma 



and of the sharp perfectly reflecting boundary confining the p lasma. While rela.xing th e 
former assumption does not change the results qualitatively (|Tvshetskiv et a/.l 12012(a ). 
relaxing the latter assumption may change the results dramatically. Firstly, the smooth 
boundary leads to a new resonant damping of surface oscillations, dramatically increasing 
the exponential damping rate jTsI in p.ip ( Marklund et a/.ll2008() . Secondly, allowing for 
boundary smoothness (with a simultaneous account for the quantum tunneling, as they 
both have the same spatial scales) should change the analytic properties of p(f2, K) in 
the lower semiplane Im(r2) < of complex fl plane, and thus may change its branch cuts 
and their contribution into (12. 7|) . This, however, is a more elaborate problem, which is 
beyond the scope of this paper and is left for future work. 



Appendix A. Branching points of ({Q,K\\) 

Let us show that the function ({il,K^^), analytically continued into Im(il) ^ 0, has 
branching points at 51 = ±51^ (Xn) S M, with the corresponding branch cuts going down 
from these two points. We start from C(r2, i^n) defined by Eq. (I2.14p for Im(J7) > 0, and 
then continuously change Im(J7) to negative values. In this process, we must consider 
how the singularities of t he function [ir^gffi,i(r)] ^ under the integral in ()2.14p change 



in the complex Kx plane (jTvshetskiv et aLll20126l ). These singularities are: 



(a) Branch cuts of the complex square root \/ ~^ defined by two parametric 
equations: 



Kx = ±i\ Kl + T, with K» > 0, r e [0, +oo). (A 1) 



(&) Branch cut of the complex logarithm in ()2.15p . taken along the negative real axis 
of the argument (SI + \/?,K)/{ri — \/?iK). This branch cut maps into two branch cuts of 
[if^ e(Sl, K)\^ in the complex plane, given by two parametric equations: 




Kx=±i\l—[- — -] - with e C, i^ii > 0, T e [0,+oo). (A 2) 

(c) Two poles Kx = {K\\ > 0) at the roots of = 0, lying symmetrically above 
and below the real axis of the complex Kx plane. 

(d) Two poles ±K^ e C at the roots of s{n,K) = 0. Note that for any K e R, 
£(0, -RT) = does not have roots with Im(i7) > 0, if the plasma equilibrium is sta- 
ble (Penrose 1960l) . which is the case considered here; therefore, for any Im(fi) > the 



poles ±KJ^ are located away from the real axis of the complex Kx plane, and thus do 
not lie on the integration contour in (|2.14l) . 

In the process of analytic continuation, as Im(r2) — > 0+, these singularities deform/move 
in the complex Kx plane, as shown in Fig. |31 We have the following cases: (i) |Re(r2)| < 
|^^„(Js^l|)|, and (ii) |Re(S2)| > |rj^(i4r||)|, where ±fly{K\\) e R are defined by the equa- 
tion p.Sp . The difference between these two cases is that for |Re(r2)| > |fii,(_fir||)|, the 
poles ±K^ cross the real axis in Kx plane and deform the integration contour when 
Im(17) 0+ and beyond to negative values, while for |Re(52)| < |fl„(Jsr||)| they do not. 
Thus we have that in these two cases, the integration contours in the function (^(D,,K\\) 
continued to Im(S7) < are different, and thus the values of (^{rt,K\\) in the lower semi- 
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Figure 3. (Color online) Singularities of the function [K'^ e(n, K)] ^ in the complex Kx plane, 
and their modification from lm(f2) > (left panels) to the limit Im(r2) — > 0-1- (right panels), 
for < Re(r2) < Qy{K\^) (upper row) and Re{Q) > fly{K^\) > (lower row), with Qy{K\\) 
defined from ([33}. The branch cuts (|A 1|) and ([A 2[l are shown with the black dashed lines and 
the blue dot-dashed lines, respectively. The poles where e{fl, K) = are shown with the 

filled circles. The arrows show the direction of motion of the singularities when Im(f2) — > 0+. 
The contour of Kx integration in (I2.14p is shown with the solid red line. 



plane of complex fl are also different for |Re(J7)| < |f2i,(iir||)| and for |Re(f2)| > |r2i,(_fir||)|. 
Hence, the points il = ztilv{K^\) e M separating these two cases must necessarily be the 
branching points of C(^^, ) in Im(f2) ^ 0, with branch cuts (separating the different val- 
ues of the analytically continued C) going down into the lower semiplane of com plex Q. At 
the br anching points zLfly{K^\), the function C(r2, K^\) has a singularity f|Tvshetskiv et al. 
20126h . 



Appendix B. Approximation of /5(f],K) in the upper vicinity of the 
branching point il. = +ri^(iir ) 

The function p(f2,K) defined in (|2.9|) contains integrals of the form 
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which need to be approximated in the upper vicinity of the branching point — Qy. The 
main contribution into the integrals (|B 1|) is from the vicinity of Kx = 0. The expansion 
of e( fl, K) near J7 = f2„, — under the integral is 

e{n,K) = {n - n^) l3{Kn) + a{Kii)Kl + O \{n - n^f] +0{K^) , (B2) 

where 

13{K\\)^ e'^{^,K)\^^^, ,,^^^, (B4) 

m) = eni^,K)\^=n., (B5) 

and the primes denoting partial derivatives with respect to the corresponding variables, 
e.g., = de/dfl. Here we have taken into account that e{i^v,K\^) — (by definition of 
and e'j,^in,K)\^^^^ = [{KjK)e'^]\^^^, = 0. 
The function C(^i in the upper vicinity of = ily(K\\) is then 

1 1 /■+°° dKx 



C(r!,i^ii 



1 1 1 



(B( 



2 27^11 ^a(if||)/3(/^||) Vf^-r!. 

Here we neglected 0[(f2 — ilt,)^] (as we are considering a small vicinity of 57 — 
0{K^) (as the main contribution into the integral is from ~ 0), and 0[{il — fty)K'^] 
(due to the combination of the above two reasons). Similarly, in the upper vicinity of 
fl = riy^K^^) we obtain 

r^dK' ^^"'^^^ - w(».,K)i^^,„ 

J-oo ''K'h{n,K') ^a(i^||)/3(if||)i^2^jjT7Ti;' 

where we have assumed that the function /(fi, K) varies slowly near the point il = f2„, 
and does not have branch cuts. 

Combining (jB 6|) and (|B 7p in (|2.9p . we obtain the approximation (|3.6I) for p(f7,K) in 
the upper vicinity of — Qy. 
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